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ABSTRACT 


Y 

The  following  cype  of  problem  arises  in  practice:  in  a  node-weighted 
graph  G,  find  a  minimum  weight  node  set  that  satisfies  certain  conditions 
and,  in  addition,  induces  a  perfectly  matchable  subgraph  of  G.  This  has  led 
us  to  study  the  convex  hull  of  incidence  vectors  of  node  sets  that  induce 
perfectly  matchable  subgraphs  of  a  graph  G,  which  we  call  the  perfectly 
matchable  subgraph  polytope  of  G.  For  the  case  when  G  is  bipartite,  we 
give  a  linear  characterization  of  this  polytope,  i.e.,  specify  a  system  of 
linear  inequalities  whose  basic  solutions  are  the  incidence  vectors  of 
perfectly  matchable  node  sets  of  G.  We  derive  this  result  by  three  different 
approaches,  using  linear  programming  duality,  projection,  and  lattice 
polyhedra,  respectively.  The  projection  approach  is  used  here  for  the  first 
time  as  a  proof  method  in  polyhedral  combinatorics ,  and  seems  to  have  many 
similar  applications  Finally,  we  completely  characterize  the  facets  of 
our  polytope,  i.e.,  ^*F  separate  the  essential  Inequalities  of  our  linear 
defining  system  from  the  redundant  ones.  ^ — 


THE  PERFECTLY  MATCHABLE  SUBGRAPH  POLYTOPE 


OF  A  BIPARTITE  GRAPH 
by 

Egon  Balas  and  William  Pulleyblank 

1.  Introduction 

Given  a  graph  G  ■  (V,E) ,  it  is  often  of  interest  to  identify  those 
node  sets  of  G  that  are  perfectly  matchable,  i.e.,  those  S£V  such  that 
<  S  >  ,  the  subgraph  of  G  Induced  by  S,  has  a  perfect  matching.  We  call 
the  convex  hull  of  the  incidence  vectors  of  perfectly  matchable  node  sets 
of  a  graph  G,  the  perfectly  matchable  subgraph  polytope  (IMS  polytope)  of  G. 

The  identification  of  the  perfectly  matchable  node  sets  of  a  graph  G 
would  of  course  become  much  easier  if  the  IMS  polytope  of  G  could  be  linearly 
described,  i.e.,  if  one  had  a  system  of  linear  inequalities  whose  basic  solu¬ 
tions  are  precisely  the  extreme  points  of  the  IMS  polytope  of  G.  The  existence 
of  such  a  linear  system  follows  from  the  by  now  classical  result  that  the 
convex  hull  of  a  finite  set  of  points  in  Kn  is  the  intersection  of  a  finite 
number  of  halfspaces  ia Bn,  i.e.,  the  solution  Bet  of  a  finite  system  of 
linear  inequalities  in  n  variables.  But  the  identification  of  such  a  linear 
system  defining  a  polytope  given  by  the  set  of  its  extreme  points  (that  are 
either  explicitly  listed  or  specified  by  some  definition,  like  here)  is 
usually  a  hard  task,  which  has  so  far  been  solved  only  for  a  few  cases.  In 
this  paper  we  give  such  a  linear  characterization  of  the  MS  polytope  of 
a  bipartite  graph.  The  case  of  a  general  graph  will  be  addressed  in  another 
paper. 

The  question  examined  here  arose  in  the  context  of  a  real  world 
problem  that  had  to  do  with  the  optimal  scheduling  of  drivers  for  a  municipal 


bus  company.  This  particular  application,  which  gave  the  initial  motiva¬ 
tion  for  our  research,  is  described  in  section  2  of  the  paper.  Section  3 
introduces  the  system  of  linear  inequalities  defining  the  IWS  polytope  of 
a  bipartite  graph  and  gives  a  first  proof  of  the  validity  of  this  linear 
characterization,  based  on  linear  programming  duality  theory.  Section  4 
gives  an  alternative  proof,  using  a  projection  technique  that  is  of  interest 
in  itself,  since  it  may  serve  as  a  proof  method  in  situations  analogous  to, 
but  different  from,  the  one  examined  here.  Finally,  section  5  gives  a  third 
proof,  based  on  the  theory  of  lattice  polyhedra. 

Section  6  of  the  paper  focuses  on  the  question  of  redundancy  in  the 
system  introduced  in  section  3,  and  gives  a  complete  characterization  of 
the  facets  of  the  FMS  polytope  of  a  bipartite  graph. 

2.  Motivation:  A  Bus  Driver  Scheduling  Problem 

The  following  problem  was  brought  to  our  attention  by  Mr.  A.  Foe6 
of  the  Operations  Research  group  of  Nederlandse  Spoorwegen,  the  Dutch 
Railway  Company. 

A  municipal  bus  company  had  to  schedule  the  tours  of  duty  of  its 
drivers,  so  as  to  cover  a  daily  set  of  trips  to  be  executed.  A  set  covering 
approach  was  used,  i.e.,  the  problem  was  formulated  as 

min{cx|Ax  >  e,  x  e  £0, 1 }n) , 

where  A  is  an  m  X  n  0-1  matrix  whose  jC^  column  represents  a  potential 
daily  (tour  of)  duty  for  a  driver,  with  a^  =  1  if  duty  j  covers  trip  i, 

«  0  otherwise,  while  c^  is  the  cost  of  duty  j,  and  e  =  (1,...,1).  In 
a  typical  case  the  matrix  A  had  about  150-200  rows  and  3000-4000  columns. 

However,  the  way  the  columns  of  A,  i.e.,  the  potential  duties,  were 
generated,  suggested  another  approach.  Initially,  a  set  of  "early  parts" 
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(morning  half-tours)  and  "late  parts"  (afternoon  half-tours)  of  duty  were  generated 
independently  of  each  other,  then  all  the  compatible  early  part-late  part 
pairs  were  explicitly  generated  as  potential  full  day  duties.  The  number 
of  early  parts  and  late  parts  was  typically  about  150  and  200  respectively, 
and  the  3-4000  columns  of  A  arose  from  the  fact  that  only  10-13%  of  the 
30,000  pairs  were  compatible  (because  of  starting  and  ending  properties  in 
space  and  time).  If  the  number  of  early  parts  and  late  parts  is  n^  and  n2, 
respectively,  and  the  ratio  of  compatible  early  part-late  part  pairs  to  all 
such  pairs  is  r,  then  n  *  r  X  n^  X  n^;  i.e.,  n  is  usually  much  larger  than 
^  +  n2- 

11  2  2 

Now  let  A  *  ^  “  (a^,  j)  b®  m  x  n^  and  m  x  n2  matrices, 

respectively,  defined  by 

1  if  early  part  j  covers  trip  i 

0  otherwise 


and 


2 

“ij 


1  if  late  part  j  covers  trip  i 

„  0  otherwise, 


1  2 

and  let  c  and  c  be  the  cost  vectors  of  early  parts  and  late  parts, 
respectively.  Further,  let  G  *  (V^UVj,  E)  be  the  bipartite  graph  whose 
node  sets  and  V2  correspond  to  the  early  parts  and  the  late  parts, 
respectively,  and  whose  edges  correspond  to  compatible  early  part-late 
part  pairs.  Then  the  above  problem  can  be  reformulated  as  follows: 

1  2 

(i)  Find  x  e{0,l)  and  x  c(0,l}  to 
(2.1)  minimize  cSc*  +  c2x2 


subject  to 

(2.2)  A^1  +  A2x2  >  e 
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and 

1  2 

(2.3)  (x  ,  x  )  is  the  incidence  vector  of  some  SEV^L'l^  such  that 
<  S  >  has  a  perfect  matching. 

(ii)  Find  a  minimum-weight  perfect  matching  in  the  graph  <  S  > 
with  edge-weights 

Sij*ci  +  cJ*  a,j)eE. 

Here,  as  before,  <  S  >  denotes  the  subgraph  of  G  induced  by 
the  node  sec  S. 

Problem  (ii)  is  of  course  polynomially  solvable;  whereas  problem  (i) 
replaces  the  original  3-4000  variable  set  covering  problem  by  a  350-variable 
see  covering  problem  with  side  condition  (2.3). 

The  solvability  of  the  problem  thus  hinges  on  whether  one  can  con¬ 
veniently  represent  condition  (2.3). 

3.  A  Linear  Characterization  of  the  FMS  Polvtope 

Let  G  *  (V^UV2,  E)  be  a  bipartite  graph  with  parts  and  V2,  i.e., 
with  node  set  V  *  an<i  e<*8e  8et  E  suc^  c*iat  every  eeE  Joins  some  node 

of  Vj  to  some  node  of  V^. 

Let  fll(G)  be  the  family  of  perfectly  matchable  node  sets  of  G,  i.e., 
771(G)  *  {Scv|<  S  >  has  a  perfect  matching}. 

For  any  SCV,  the  incidence  vector  (characteristic  vector)  of  S  is 
x  e[0,l}^  such  that  x^  *  1,  JcS,  x^  ■  0,  jcV\S.  Let  2(G)  be  the  set  of 
incidence  vectors  of  members  of  7a(G),  and  for  any  set  T,  let  conv  T  denote 
Che  convex  hull  of  T. 

Our  objective  in  this  section  is  to  give  a  linear  system  of  inequali¬ 
ties  defining  conv  %(G) ,  i.e.,  the  WS  polytope  of  G. 
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Whenever  it  is  not  confusing,  we  will  write  V{  for  ffj(G)  and  Z  for  Z(G) . 
Many  problems  Involving  matchings,  in  particular  in  bipartite  graphs, 
can  be  shown  to  be  special  cases  of  certain  matroid  problems.  For  instance, 

if  G  *  ^  C^e  bipartite  graph  introduced  above  and  for  k  =  1,  2, 

is  the  family  of  those  edge  sets  that  meet  every  node  in  at  most 

once,  then  the  system  “  (E,  is  a  matroid;  and  the  intersection  of 
the  two  matroids  and  Mj  is  the  independence  system  (E,  ^f)*^),  v^lere 
is  simply  the  family  of  all  (not  necessarily  perfect)  matchings 
in  G.  The  matching  polytope  of  G  is  then  the  convex  hull  of  incidence 
vectors  of  all  members  of 

Another  example,  more  closely  related  to  our  problem  is  the  following. 
In  an  arbitrary  graph  H  with  node  set  N,  let  J  be  the  family  of  those  sub¬ 
sets  of  N  covered  by  some  matching.  Then  the  system  (N,  J) ,  as  shown  by 
Edmonds  and  Fulkerson  [4] ,  is  a  matroid. 

In  such  cases  as  the  above,  results  on  matroid  polyhedra  due  to 
Edmonds  [2,  3]  lead  to  linear  characterizations  of  the  type  that  we  are 
interested  in.  However,  these  results  are  not  applicable  to  our  case,  since 
the  IWS  polyhedron  of  a  graph  (bipartite  or  not)  does  not  have  a  matroidal 
structure.  To  see  this,  it  is  sufficient  to  recall  the  fact  that  every 


S  c 771(G)  is  of  even  cardinality. 

We  now  briefly  state  our  notational  conventions.  An  edge  joining 
nodes  i  and  j  is  denoted  (i,j).  For  S,  TCV,  the  set  of  edges  joining  nodes 
in  S  to  nodes  in  T  is  denoted  (S,T).  For  SCV,  I"(S)  denotes  the  set  of 


nodes  adjacent  to  some  node  in  S.  Clearly,  if  SCV^,  then  and 

vice  versa.  For  the  sake  of  brevity,  we  write  T(i)  for  F({i]). 

For  any  x  and  any  SCV,  we  let  x(S)  =  Kx^ieS). 
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Next  we  state  the  linear  system  defining  the  WiS  polytope  of  G,  i.e., 
the  convex  hull  of  X. 

Theorem  3.1.  Let  P  be  the  convex  polytope  consisting  of  those 

Ivl 

x  e  1R  '  satisfying 

(3.1)  0  <  \  <  1.  icV 

(3.2)  x(Vx)  -  x(V2)  -  0 
and 

(3.3)  x(S)  -  x(r(s))  <  o,  *scvr 

Then  P  *  conv  X. 

Proof.  It  is  easy  to  see  that  conv  XCP.  For  let  x  be  any  vertex 
of  conv  X;  then  x  is  the  incidence  vector  of  some  T  e#!,  hence  (3.1)  holds 
trivially.  Further,  (3.2)  is  the  requirement  that  ItOV^  *  \tHV2|,  and 

(3.3)  simply  states  that  for  any  SCV^,  T  must  contain  at  least  as  many 
nodes  of  T(S)  as  of  S.  Both  of  these  requirements  are  readily  seen  to  be 
necessary  conditions  for  <  T  >  to  have  a  perfect  matching,  and  together  they 
constitute  the  "easy"  part  of  the  well-known  KBnig-Hall  theorem  [11],  [6]. 

To  prove  the  converse,  namely  that  P £ conv  X,  we  will  show  that 
every  vertex  of  P  belongs  to  X.  This  will  be  done  by  showing  that  for 
any  vector  c  »  (c^zieV)  of  real  node  costs  there  is  an  optimal  solution  x* 
to  the  linear  program 

(L)  max{cx|x  e  P), 

such  that  x*  c  X.  Since  every  vertex  of  P  is  the  unique  optimal  solution 
to  such  a  linear  program  for  some  c,  this  will  give  the  result. 

We  define  a  vector  c  ■  (c^ : (i, J) eE)  of  edge  costs  by  letting  c  =  ^ 
for  all  (i,j)cE.  For  any  matching  MCE,  if  S  is  the  set  of  nodes  covered 
by  M,  then  M  is  a  perfect  matching  in  <  S  >,  and 
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(3.4,  I(cij:(i,j)cM)  -  K^rieS). 


Conversely,  for  any  S  efll  and  any  perfect  matching  M  in  <  S  >,  M  is  also 
a  matching  in  G,  and  (3.4)  holds.  Therefore  the  problem  of  maximizing  cx 
over  x  e  X  can  be  solved  by  finding  a  maximum-weight  matching  (in  terms  of 
the  edge-weights  c)  in  G. 

Let  M*  be  such  a  matching,  and  let  x*  be  the  incidence  vector  of  the 
node  set  S*  covered  by  M*.  We  will  show  that  x*  is  an  optimal  solution  to 
the  linear  program  (L) ,  by  constructing  a  feasible  solution  to  the  dual  of 
(L)  having  the  same  objective  function  value  as  (L) . 

Since  edge-variables  are  two-indexed,  we  amend  our  notational  con¬ 
ventions  by  writing,  for  S,  TCV,  u(S,T)  *  £(Uj.:ieS,  jeT),  and  u(i,T)  =  u({i},  T), 
u(S,  j)  =  u(S,(j}). 

The  graph  G  being  bipartite,  the  Incidence  vector  u*  of  the  matching 
M*  is  an  optimal  solution  to  the  linear  program 


max  cu 


<v 

u(i,V2)  <  1 

icVL 

u(VL,j)  <  1 

j  «V2 

u  >  0 

whose 

dual  is 

min  t(Vj)  +  t(V2) 

(Di> 

h  +  Cj  -  Cij 

t  >  0 


Let  t*  be  an  optimal  solution  to  (D^) .  By  linear  progranming 


duality, 
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(3.5)  E(c  *  t*^)  +  t*(V2>. 

We  now  write  down  the  linear  program  (D) ,  dual  to  (L): 
min  yO^)  +  y(V2> 


(3.6) 

yi 

+  E(*s:SfiV1# 

ieS)  >  Cj 

ieV 

(3.7) 

yj 

-  E<VS=V 

3*  r(s))  >  Cj 

JCV; 

(3.8) 

yi 

»  Tj  >  0» 

i*V^,  jeV2 

(3.9) 

ZS 

>  0,  s  £  Vx 

;  zy  unconstrained. 

Now  let  y*  •  tj  for  ieV^,  y^  =  t *  for  jeV2;  and  zg  *  0  for  all  ScV^. 
Then  (3.8)  -  (3.9)  are  satisfied,  and 

(3.10)  y*(Vx)  +  y*(V2)  =  t*(vp  +  t*(V2) 

*  :  (i,j)«K*)  =  cx*. 

Next  we  will  describe  a  procedure  for  redefining  the  value  of  z  for 

u 

certain  subsets  SCV^  in  such  a  way  as  to  satisfy  (3.6)-(3.7),  without  changing 
the  value  of  any  y*.  Therefore,  the  vector  (y*p  z)  obtained  in  this  way  will 
be  the  optimal  solution  to  (D)  required  for  the  completion  of  our  proof. 

At  all  stages  of  the  procedure,  the  vector  (y*,  z)  will  satisfy  the 
following  two  synmetrlc  properties: 

(3.11)  If  for  some  i«V^ 

y*  +  I(zs:SCV1,  ieS)  =  c£  -  e 


for  some  e  >  0, 
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then  for  every  jer(i), 

y*  -  E(z  • SCV.,  je? (S)>  >  c.  +  e. 
j  a  i  J 

(3.12)  If  for  some  jeV2 

y*  -  E(z.:S5V,,  jcT(S))  =  c  -  c  for  some  e  >  0, 

3  b  1  j 

then  for  every  i«r(j), 

y*  +  E(zg:S  cvx,  ieS)  >  ^  +  e. 

These  conditions  state  that  if  the  current  solution  violates  the 
inequality  associated  with  some  node  by  an  amount  e,  there  is  a  surplus 
of  at  least  c  at  every  adjacent  node.  By  the  initial  definition  of  ( y *,  z) 
and  in  view  of  the  inequalities  t*  +  t*  >  c^ ^ ,  conditions  (3. 11) - (3. 12) 
are  satisfied  initially. 

Define 

S0  =  tieVilyi  +  r<«s:scv1#  its)  <  Cl}, 

To  =  U«V2|y*  -  E(zs:ScV1,  jef(S))  <  Cj}. 

Note  that  by  (3.11)  and  (3.12),  no  icSQ  and  jcTQ  are  adjacent. 

If  at  any  stage  of  the  procedure  S  *  T  *  0,  then  (3.6)  and  (3.7) 

o  o 

are  satisfied  and  we  are  done.  If  Sq  #  0,  let  s  =  0  and  perform  Reduction  1. 

If  S  =0  but  T  f  0,  let  t  =  0  and  perform  Reduction  2. 
o  o 

Reduction  1.  Let 

e  =  min  c,  -  y*  -  E(zc:S£V  ,  ieS) 

1151 
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and  define  zc  *  e(>0).  Then  (3.11)  and  (3.12)  are  still  satisfied 
s 

(since  (3.11)  was  satisfied  before),  but  the  set 
Ss+1  S  ^ieVi lyj  +  S(zs:S  =Vr  ieS)  < 
is  a  proper  subset  of  Sg. 

If  Sg+1  =  0,  Reduction  1  is  complete;  otherwise  set  s  -  s  +  1  and 
repeat  Reduction  1. 

Reduction  2.  Let 

c  =  min  c  -  y*  +  S(zQ:S  CV. ,  jer(S)) . 

J*Tt  j  j  SI 

Then  e  >  0.  Define  z_  »  z„  -  e,  S  *  V  \F(T  ),  and  arr  *  zr  +  e.  Note 

11  tit  &t  bt 

that  the  effect  of  this  change  is  to  decrease  c  -  y*  +  E(z  :ScV  ,  jeF(S)) 

3  3  b  1 

by  e  for  jeTt  and  to  leave  it  unchanged  for  and  also  to  decrease 

y*  +  I(ZgtSCV^,  ieS)  -  by  l  for  ic  (Tt)  but  to  leave  it  unchanged  for 
ieVj\r(Tt) . 

Conditions  (3.11)  and  (3.12)  still  hold  (since  (3.12)  was  satisfied 
before),  and  the  new  zg  still  satisfy  (3.9);  but  the  set 

Tt+1  2  "  E<VSCV1*  j«r(S))  <  Cj) 

is  a  proper  subset  of  Tt< 

If  Tt+^  *  0.  Reduction  2  iB  complete;  otherwise  set  t  •-  t  +  1  and 
repeat  Reduction  2. 

After  at  most  |SQ!  <  iterations  of  Reduction  1  and  at  most 

l^0l  £  |Vj|  Iterations  of  Reduction  2  we  obtain  a  vector  (y*,  z)  satisfying 
(3.6)-(3.9),  and  thus  the  proof  of  the  theorem  is  complete.!] 

At  this  point  some  remarks  are  in  order. 
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First,  there  is  a  certain  lack  of  symmetry  in  the  linear  system 
(3.1)-(3.3)  defining  conv  X,  in  that  it  contains  inequalities  only  for  sub¬ 
sets  S  of  V^,  but  not  for  subsets  T  of  Vj.  The  analogous  inequalities  for 
subsets  of  would  be 

(3.13)  x(T)  -  x(F(T))  <  0,  V  TCV2. 

These  are  clearly  valid  and  could  have  been  Included  in  the  system,  but 
they  can  also  be  derived  from  (3.1) -(3. 3).  For  if  TCV^  and  we  define 
S  e  Vj\r(T),  then  F(S)  CVjVT;  and  by  subtracting  (3.2)  from  the  inequality 
x(S)  -  x(T(S) )  <  0,  we  obtain  x(V2\T(S))  -  x(F(T))  <  0.  But  since 
F(S)  CV2\T  implies  TCV2\T(S),  and  since  x  >  0,  this  last  inequality 
implies  x(T)  -  x(r(T))  <  0. 

If  we  had  included  the  inequalities  (3.13)  in  our  system  defining 

conv  X,  then  Reductions  1  and  2  could  have  been  made  completely  symnetric 

by  using  the  new  dual  variables  that  would  have  been  introduced. 

Second,  suppose  is  such  that  the  graph  <  S  U  F(S)  >  is  dis¬ 

connected,  with  components  <  S^Ur(Sk)  >,  k  =  l,...,q.  Then  the  inequality 
x(S)  -  x(r(S))  <  0  is  the  sum  of  the  q  inequalities  x(Sfe)  -  x(r(Sfe))  <  0, 
k  =  l,...,q,  hence  redundant.  Now  suppose  <  SUT(S)  >  is  connected  and  K  is  the 
node  set  of  the  component  of  G  containing  <  SUT(S)  >,  with  *  KOv^,  i  =  1,2, 

but  the  graph  <  (K^\S)  U  (K2\T(S))  >  is  disconnected,  with  components  <  Tk  >, 
k  k 

k  *  1, Let  1^  =  1  HV^,  i  =  1,2.  Then  for  k  =  l,...,q,  we  have 

f(Tj  U  S)  c  T2  ur(S) ,  or  else  removing  the  node  set  SUT(S)  from  G  would  not 

make  <  T^  U  T2  >  a  maximal  connected  subgraph.  Also,  F(T^US)  ST2UF(S),  or 

else  <  T^UT2  >  would  not  be  connected.  Thus  we  conclude  that  F(T^US)  =  T^  UT(S) . 
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But  then  adding  the  q  inequalities  x(T^  US)  -x(F(Tj  US))  <  0,  k  =  l,...,q, 
and  subtracting  (q  -  1)  times  the  equation  (3.2),  yields  the  inequality 
x(S)  -  x(r(S))  <  0,  which  is  therefore  redundant. 

We  have  thus  shown  than  Theorem  3.1  remains  true  if  (3.3)  is  replaced 
by 

x(S)  -  x(r(S))  <  0  for  all  SCV^  such  that  the  graphs 
(3.3')  <  SUr(S)  >  and  C^XS)  U  (K2\T(S)>  are  connected,  where  <  K  >  is 

the  component  of  G  containing  <  SUF(S)  >  ,  and  i  =  1,2. 

Third,  note  that  if  c  is  integer  valued,  then  so  is  c”,  and  thus  t* 
can  be  chosen  to  be  integer  valued.  Then  each  iteration  of  Reduction  1  or 
2  will  result  in  integer  e  and  hence  in  integer  valued  (y*,  z) .  Thus  for 
any  integer  valued  c,  the  linear  program  (D),  dual  to  (L) ,  has  integer 
optimal  solutions.  Thus  our  linear  system  defining  the  W1S  polytope  of 
a  bipartite  graph  is  totally  dual  integral.  (This  concept  was  introduced 
by  Hoffman  [9]  and  used  extensively  by  Edmonds  and  Giles  [5] .  See  also 
Schrijver  [12].) 

Fourth,  if  we  set  c^  =  1  for  all  ieV^  and  c^  s  0  for  all  jcV2,  then 
the  value  of  (an  optimal  solution  to)  (L),  and  hence  of  (D) ,  is  the  cardi¬ 
nality  of  a  maximum  matching  in  G.  Now  suppose  G  has  no  matching  that  covers 
all  ieV;  then  if  (y*,  z*)  is  an  optimal  integer  solution  of  (D) , 

y*(Vx)  +  y*(V2)(«  max{cx|x  c  P})  <  \vx | . 

Since  each  y*  is  a  nonnegative  integer,  this  implies  that  y*  *  0 
for  some  ieV^.  But  since  (y*,  z*)  must  satisfy  (3.6),  there  must  be  some 
SCV^  such  that  z*  >  0.  Now  suppose  the  optimal  solution  (y*,  z*)  is 


13 


chosen  such  that  the  number  of  positive  components  of  2*  is  minimum,  and 
let  be  such  that  2*  >  0.  Then  |s|  >  |r(S)|;  for  if  not,  then 

by  adding  z*  to  y*  for  ieS,  subtracting  z*  from  y*  for  jcF(S),  and  then 

Si  o  j 

setting  Zg  =  0,  we  could  obtain  a  new  optimal  solution  to  (D)  with  fewer 
positive  components  of  z,  a  contradiction.  Thus  we  obtain  the  hard  part  of 
the  KBnig-Hall  Theorem,  namely  that  if  G  =  (V1  UV2,  E)  has  no  matching  that 
covers  all  of  Vp  then  there  exists  ScV^  such  that  >  |  F(S)  | .  Furthermore, 
this  last  result  combined  with  our  second  remark  gives  a  strengthened  version 
of  the  hard  part  of  the  KBnig-Hall  Theorem:  for  G  such  that  |v^[  =  |v2 |  to 
have  a  perfect  matching,  it  is  sufficient  that  the  condition  |s|  <  |r(S) |  be 
satisfied  for  every  S  CV^  such  that  <  S  UT(S)  >  and  <(K1\S)  U  (K2\p(S))>  are 
connected,  where  K  is  the  node  set  of  the  component  of  G  containing  <  SUF(S)  >, 
and  Kt  «  KflVp  i  =  1,2. 

Fifth,  any  optimal  solution  (y*,  z*)  to  (D)  can  be  seen  to  have  the 

following  property.  There  exists  a  nested  sequence  of  sets  0  #  U  CD  , C . . . 

n  n- 1 

CU^CUgCVp  such  that  for  any  SCVp  z*  >  0  if  and  only  S  =  Uj  for  some 
ic{0,.,.,n).  This  is  so  because  if  we  did  s  Iterations  of  Reduction  1, 
we  will  have  defined  sets  0  #  S  CS,  C...CS.CS  .  If  we  did  t  iterations 

S  8-1  1  O 

of  Reduction  2,  we  will  have  defined  sets  S  CS,  £...  CS  .  Further,  from 

o  l  t 

(3.11)  and  (3.12),  SqCSo.  Combining  these  sequences  gives  the  claimed 
sequence  (U^:i  *  0,1,..., n). 

Finally,  we  have  shown  that  for  any  optimal  solution  t*  to  the 
node  covering  problem  (D^),  there  is  an  optimal  solution  (y*,  z*)  to  (D) 
for  which  y*  *  t*.  Of  course  the  converse  is  also  true:  if  (y*,  z*)  is 
an  optimal  solution  to  (D) ,  then  setting  t*  =  y*  gives  an  optimal  solu¬ 
tion  to  the  node  covering  problem  (D^). 
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4.  An  Alternative  Derivation  via  Projection 
In  this  section  we  give  an  alternative  derivation  of  the  linear  system 
defining  the  R4S  polyhedron  of  a  bipartite  graph,  based  on  a  polyhedral 
interpretation  of  Benders's  partitioning  theorem  [1],  This  approach  is  of 
more  general  interest  than  its  particular  use  in  this  paper,  since  it  provides 
a  technique  for  projecting  a  polyhedron  in  En,  or  some  (not  necessarily 
polyhedral)  subset  of  a  polyhedron  in]Rn,  into  some  specified  subspace  of  ]Rn. 
To  be  specific,  let  Q  be  an  arbitrary  subset  of 1R^,  and  let 

Z  =  t(u,x)  e  lRp+<i|Au  +  Bx  <  b,  u  >  0,  x  c  Q} 

where  A,  B  and  b  are  m  X  p,  m  X  q,  and  m  X  1  matrices,  respectively,  such 
that  Z  f  0.  The  pro lection  of  Z  into  the  subspace  of  the  x-variables  is 
defined  as 

X  *  {x  elR^I  there  exists  u*Rp  such  that  (u,x)cZ}. 

We  are  interested  in  describing  the  set  X  in  a  way  similar  to  Z,  l.e.,  by 

a  set  of  linear  inequalities  plus,  of  course,  the  condition  x  c  Q.  The 

following  theorem  accomplishes  this. 

Before  stating  the  result,  we  recall  that  a  polyhedral  cone  C  is  the 

intersection  of  a  finite  number  of  halfspaces  through  the  origin,  and  a 

pointed  cone  is  one  of  which  the  origin  is  an  extreme  point.  A  ray  of  a 

cone  C  is  the  set  R(y)  of  all  nonnegative  multiples  of  same  ycC,  called 

the  direction  (vector)  of  R(y).  A  vector  y«C  is  extreme,  if  for  any  y^, 
112  12 

y2eC,  y  ■  y(y  +  y  )  implies  y  ,  y  *R(y) .  A  ray  R(y)  is  extreme  if  its 
direction  vector  y  is  extreme.  A  pointed  polyhedral  cone  has  a  finite 
number  of  extreme  rays,  and  is  the  conical  hull  of  its  extreme  rays. 

Of  course,  for  every  nonsero  x  c  R(y),  we  have  R(x)  -  R(y)  and  consequently 
every  cone  that  contains  more  than  the  origin  has  an  infinite  number  of 
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extreme  direction  vectors.  However  the  smallest  set  of  vectors  of  which 
a  cone  is  the  conical  hull,  consists  of  one  direction  vector  from  each 
extreme  ray. 

For  a  cone  C  we  let  extr  C  denote  such  a  (finite)  set  of  extreme 
direction  vectors.  Note  that  £?tr  C  is  uniquely  determined  up  to  positive 
multiples. 

Theorem  4.1.  Let  Z  and  X  be  defined  as  above,  and  let 
W  *  (ve  ]Rm(vA  >0,  v  >  0}  . 

Then 

X  *  {x  c!Rqj(vB)x  <  vb,  ¥  ve  extr  W;  x  e  Q}. 

Proof.  The  polyhedral  cone  W  is  a  subset  of  hence  pointed. 
Therefore  W  is  the  conical  hull  of  its  extreme  rays,  and  any  x  s  ]Rq  satisfies 
the  inequality  (vB)x  <  vb  for  every  extreme  direction  v  of  W,  i*  and  only 
if  it  satisfies  it  for  all  veW. 

Now  let  x  e  X;  then  x  c  Q  and  there  exists  uelp  such  that  u  >  0 
and  Au  +  Bx  <  b.  Further,  let  veW;  then  vBx  <  vb  -  vAu  <  vb,  since  u  >  0  and 
vA  >  0  imply  vAu  >  0.  Thus  (vB)x  <  vb,  ¥  vc  extr  W. 

Conversely,  suppose  x  e  ]Rq  satisfies  x  e  Q  and  (vB)x  <  vb,  ¥  vc  extr  W. 
Then  there  exists  no  ve  ]Rm  such  that  vA  >  0,  v  >  0  and  v(b  -  Bx)  <  0. 
Therefore,  from  Farkas's  well  known  Letmna,  there  exists  some  ue  ]RP 
such  that  u  >  0  and  Au  <  b  -  Bx.  But  then  x  c  X.|j 

Note  that,  if  W  ■  {0}  (like  for  instance  in  the  case  when  A  <  0) ,  then 
X  -  (x  e  Eq|x  e  Q}. 

We  now  turn  to  our  problem  of  giving  a  linear  characterization  of 
the  ms  polytope  of  a  bipartite  graph  G.  Although  we  are  looking  for  a 
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linear  system  in  terms  of  the  variables  xi  associated  with  the  nodes  of  G, 
we  will  start  with  the  much  easier  task  of  giving  a  linear  character¬ 
ization  in  terms  of  variables  associated  with  both  nodes  and  edges.  Such 
a  linear  system  of  course  defines  a  polyhedron  in  a  higher  dimensional 
space  than  the  one  that  we  are  looking  for,  however  by  projecting  this 
polyhedron  into  the  space  of  the  node  variables  we  will  obtain  the  system 
of  Theorem  3.1. 

Recall  that  the  IMS  polytope  of  G  is  conv  Z,  where  Z  is  the  set  of 
incidence  vectors  of  perfectly  matchable  node  sets  of  G.  Let,  as  before, 
a  variable  x^  be  associated  with  node  i  of  G,  and  let  a  variable  u^  be 
associated  with  edge  (i,j)  of  G.  As  in  section  3, we  write  u(S,T)  «  Ku^iicS,  j*T), 
u(i,T)  =  u((i],T),  and  u(S,j)  -  u(S,{j}). 

Ivl 

It  is  not  hard  to  see  that  a  0-1  vector  x  eK1  1  is  the  incidence 

vector  of  some  perfectly  matchable  node  set  of  G  if  and  only  if  there  exists 
|e| 

some  integer  uclR*  1 ,  such  that 

ua.ru))  -  xt  -  o  uv1 

(4.1)  u<r<j),j>  -  x  =  o  j«v2 

uij  >  o,  (i,j)cE. 

Furthermore,  since  the  coefficient  matrix  of  (4.1)  is  totally 
unimodular,  the  integrality  condition  on  u  can  be  omitted,  and  the  0-1 
condition  on  x  can  be  replaced  by 

(4.2)  0  <  xjL  <  1,  isV. 

Thus  (4.1)  and  (4.2)  provide  a  linear  characterization  of  conv  X 
in  terms  of  node  and  edge  variables.  One  way  of  obtaining  a  linear  char¬ 
acterization  in  terms  of  the  node  variables  only,  is  then  to  project 
the  polyhedron  defined  by  (4.1),  (4.2)  into  the  subspace  of  the  node  variables. 


17 


To  this  end,  we  first  rewrite  (4.1)-(4.2)  as  a  system  of  linear 
inequalities.  This  can  be  done  in  several  ways,  and  we  choose  to  (a)  change 
the  sign  of  the  equations  jeV^;  ($)  replace  all  equations  by  inequalities 
of  the  form  <;  and  (y)  add  all  the  inequalities  thereby  obtained  for  icV^ 
and  j*V2,  and  change  the  direction  of  the  resulting  inequality.  This  yields 
the  system 


-u(i,r(i)) 

+  xi 

< 

0 

ieV1 

u(r(j),j> 

'xj 

< 

0 

j  ®V2 

-x(Vx)  + 

*(v2) 

< 

0 

uij 

> 

0 

0  £  xi 

< 

1 

ieV 

which  is  equivalent  to  (4.1)-(4.2).  Note  that  the  coefficient  matrix  of 
(4.3)  is  still  totally  unimodular. 

We  now  apply  Theorem  4.1  to  this  system.  The  set  Q  and  the  matrices 
A,  B  and  b  that  define  Z  of  Theorem  4.1  are  in  this  case  as  follows: 

Q  *  {x  s!rIV1|  -  xO^)  +  x(V2)  <0,  0  <  x£  <  1,  ieV}; 

A  is  the  node-edge  incidence  matrix  of  G,  with  the  signs  of  the 
rows  indexed  by  changed; 

B  is  a  diagonal  matrix  of  order  jv|,  with  +1  for  the  diagonal  entries 
indexed  by  V^,  and  -1  for  those  Indexed  by  V2;  and,  finally, 

b  is  the  0  vector  with  |v|  components. 

Now  the  cone  W  of  Theorem  4,1  is 


W  »  < 

-vi  +  Vj  >  0, 

j*V2,  (i, j) SE 

l 

\  >  0, 

ieV 

r 


T 


1 


i 

( 


y 
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and  in  order  to  project  the  polyhedron  defined 
of  the  node  variables,  we  have  to  characterize 
Theorem  4.2.  The  vector  v*W  is  extreme 
a  >  0  such  that  either 

f O'  for  exactly  one  i  =  j^«V2 


(4.4) 


(.0  for  all  ieV^V^jJ 


by  (4.3)  into  the  subspace 
the  extreme  rays  of  W. 
if  and  only  if  there  exists 


or 


(4.5) 


its  u  r(s) 

otherwise 


for  some  SCV^  such  that  <  S  U  F(S)  >  is  connected. 

Proof.  Sufficiency.  Let  v«W  be  of  the  form  (4.4),  and  assume  for 

11  2 

the  sake  of  contradiction  that  v  is  not  extreme,  i.e.,  v  *  j(v  +  v  )  for 
some  v*,v2eW\R(v) .  Then  v*  *  v2  *  0,  ¥  itV^  U ,  and  v1,v2«R(v).  Thus 

v  is  extreme. 

11  2 

Now  let  veW  be  of  the  form  (4.5),  and  again  assume  that  v  =  — (v  +  v  ) 
for  some  v\v2eW.  Then  *  v2  *  0  for  ie(V^\S)  U  (V2\T(S)),  and 

(4.6)  v*  +  v2  «  2o,  icSUF(S). 

Note  that  from  (4.6),  for  any  i«S,  je"(S),  v*  >  v*  if  and  only  if 

2  2  k  k 

vi  <  Vjj  but  the  constraints  of  W  imply  v^  <  v^ ,  k  *  1,2,  for  any  such  pair 

v  k  -  - 

i,j.  Hence  v£  *  v^,  k  *  1,2,  for  all  pairs  i«S,  jtI"(S);  and  since  <  SU  *(S)  > 

k  k  k 

is  connected,  it  follows  that  v^  -  v^  *  a  (constant),  k  *  1,2,  for  all 

1  2 

i,jeSU  T(S).  Therefore  v  ,  v  eR(v),  i.e.,  v  is  extreme. 

Necessity.  Let  v  be  an  extreme  vector  of  W,  and  let  T  *  (icVjv^  >  Oj. 
Define  S  »  TflV^,  and  consider  first  the  case  where  S  *  (S.  Then  if 


I 

I 
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T  =  with  t  >  1,  and  if  denotes  the  unit  vector  inR^  with  1 


in  position  j,  we  have 


v  -  v.  e .  +  . . .  +  v.  e, 

h  h  ^t 


1,  1  .  2. 

*  2(v  +  v  i • 

where  v1  *  2v  e  ,  v2  «  2(v  e  +  ...  +  v.  e.  )f  with  v1,  v2sW,  and  v1iR(v) . 

22  J1  J2  J2  ^t 

v  £R  (v) .  Thus  if  j T J  >  1,  v  is  not  extreme,  contrary  to  the  assumption.  We 

conclude  that  if  S  -  0,  then  |t|  *  1  and  thus  v  is  of  the  form  (4.4). 

Now  consider  the  case  when  S  +  0.  Then  r(S)cTnV-,  or  else  there 

exists  ieS,  jcr(S)  such  that  >  0,  v  =  0,  i.e.,  v  violates  some  constraint 

of  W.  Also,  T(S)  2 T  1 V. ,  or  else  there  exists  jQeV2\r(S)  such  that  v  >  0. 

.  -  ^  o 

12 

But  then  for  any  c  satisfying  0  <  c  <  v  ,  the  vectors  v  and  v  ,  obtained 

1  ~  J°  2 

from  v  by  replacing  with  v  *  v  +  c  and  v  *  v  -  *,  respectively, 

^o  ^o  ^o  ^o  ^o 

112  12 

satisfy  the  equation  v  =  -j(v  +  v  alt^ou8h  v  »v  *M'\R(V)»  contrary  to  the 

assumption  that  v  is  extreme.  We  therefore  have  T(S)  *  TOV^,  i.e.,  T  =  S  UT(S), 

We  claim  that  <  T  >  is  connected.  For  suppose  not,  and  let  K  be  the 

11  2 

node  set  of  a  component  of  <  T  >.  Then  v  *  j(v  +  v  ) ,  where 


i«V1  UVjNjR 


2vt  IcK 


ieVj^  U  VjNjC, 


1  2 

while  at  the  same  time  v  ,v  *W\R(v),  contrary  to  the  assumption  that  v  is 
extreme.  Thus  <  T  >  ■  <  SUT(S)  >  is  connected. 

Finally,  we  claim  that  v ^  -  a,  i*T,  for  some  constant  q  >  0.  For 
suppose  not;  then  again  v  -  ^(v1  +  v2) ,  with  v1  and  v2  defined  by 


20 


and 


min{v^: jeT} 
0 


ieT 

ieV  UV2'\T 


2v^  -  min{Vj:jcT] 
0 


ieT 

ieVx  UV2'\T 


1  2 

while  v  ,v  eW\R(v),  contrary  to  the  assumption  that  v  is  extreme. 

This  proves  that  if  S  #  0,  then  v  is  of  the  form  (4.5).jj 

Having  described  the  extreme  rays  of  W,  we  can  now  apply  Theorem  4.1 
to  the  system  (4.3).  The  extreme  direction  vectors  of  the  form  (4.4)  give 
rise  to  inequalities  >  0,  i«V2»  which  are  redundant  (since  they  are  part 
of  the  definition  of  Q) .  The  extreme  vectors  of  the  form  (4.5)  give  rise 
to  an  inequality  x(S)  -  x(I*(S))  <  0  for  every  SCVj^  such  that  <  S  J  f(S)  > 
is  connected. 

If  G  is  connected,  then  the  inequality  x(Vp  -  x(V2)  <  0,  which  can 
also  be  written  as  x(Vj)  -  x(I(Vj))  <  0,  obtained  from  the  extreme  vector 
of  W  that  corresponds  to  S  =  V^,  together  with  the  inequality  -  x(V^)  +  x(V2)  <  0 
of  (4.3),  gives  rise  to  the  equation  x(V^)  -  x(V2)  =0.  If  G  is  disconnected 
with  components  <  >  ,...,  <  Kfc  >,  where  =  S^UT(S^),  i  =  l,...,t,  then 

the  equation  x(V^)  -  x(Vj)  *  0  is  obtained  by  first  adding  the  inequalities 
x(S^)  -  x(f(S^))  <  0,  i  =  l,...,t,  and  then  combining  the  resulting  inequality, 
x(V^)  -  x(V2)  <  0,  with  the  inequality  -  x(V^)  +  x(V2)  <  0  of  (4.3). 

Thus  applying  Theorem  4.1  to  the  system  (4.3),  we  obtain  the  linear 
characterization  of  the  R4S  polytope  of  G  given  in  Theorem  3.1,  except  for 
those  inequalities  (3.3)  such  that  <  SU  KS)  >  is  disconnected,  which  are 
missing.  But  these  inequalities  are  redundant,  as  shown  in  the  remarks 
following  Theorem  3.1,  where  the  system  (3.3)  was  replaced  by  (3. 3'). 
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5.  A  Third  Derivation  via  Lattice  Polyhedra 
Lattice  polyhedra  were  introduced  by  Hoffman  and  Schwartz  [10] 

(see  also  [7],  [8])  as  a  class  of  Integer  polyhedra  that  generalizes  both 
matroid  polyhedra  and  bipartite  matching  polyhedra.  We  will  show  that  the 
mS  polytope  of  a  bipartite  graph  can  also  be  expressed  in  this  form. 

A  lattice  £  is  a  partially  ordered  set  closed  under  two  associative 
and  comnutative  binary  operations,  A  and  V,  and  such  that 


(5.1) 


for  a,b«sd  ,  a  A  b  <  a,b  <  a  v  b; 


a  <  b  »  a  =  a  b,  b  *  a  V  b. 


To  define  a  lattice  polyhedron,  we  further  need  a  set  V.  and  a 

V 

mapping  -  2  that  satisfies  for  every  W3  es£, 

(5.2)  WX<W2<W3  implies  f  (W^  f!  f  (»3)  s  f  (wp 

(5.3)  f(wp  0f(W2)  Cf(Wx  VW2)flf(W1  A  w2) 

(5.4)  f(W1)Uf(W2)cf(W1  V  W2)  Uf(Wx  A  w2). 

and  a  submodular  function  r:s£ -•  TL  +  (the  set  of  nonnegative  integers).  The 

basic  result  on  lattice  polyhedra  [10]  can  then  be  stated  as  follows. 

I'l/I 

Theorem  5.1.  ?or  any  nonnegative  integer  deS1  *,  the  convex  poly- 

t"/| 

hedron  whose  points  are  those  x  sK1*'1  satisfying 

(5.5)  0  <  x  <  d 


(5.6)  Z(xi:ief(W))  <  r(W),  VWesd, 

has  only  integer  vertices.  Moreover,  the  linear  system  (5.5),  (5.6)  is 
totally  dual  Integral. 
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To  apply  this  theorem  to  our  case,  we  let  I  be  the  collection  of  all 
WCV2  ordered  by  set  inclusion,  and  we  define  the  operations  v  and  A  to 
be  U  and  0  ,  respectively.  Then  £  is  well  known  to  be  a  lattice.  We  let 
V.  =  V,  the  node  set  of  G. 

For  tftsf  we  define  f(W)  =  S  UW,  where  S  =  r_1(W)  is  the  maximal 
subset  of  Vx  such  that  T(S) =W.  Equivalently,  S  consists  of  all  those 
nodes  of  adjacent  only  to  nodes  in  W. 

Now  for  W i  es£,  i  =  1,2,3,  condition  (5.2)  requires  that  CW^ 

imply 

(W1us1)  n  (W3us3)  C  (W2us2), 

where  =  T1^),  i  *  1,2,3.  Since  implies  S1CS2CS3.  this 

condition  is  satisfied. 

Further,  for  i  *  1,2,  (5.3)  requires  that 

(w1us1)  n  (w2us2)  cwlnw2ur"l(w1nw2), 

where,  again,  S^T"1^),  i  -  1,2.  Since  (Vx  USj)  fl  (»2  U  S2)  =  (W^l^)  U  (Sx  OS  ), 
and  since  it  is  easily  checked  that  S^ftS2  =  r  ^(W^  H W2 ) ,  this  requirement 
is  also  satisfied. 

Finally,  for  Wt  «£,  i  «  1,2,  (5.4)  requires  that 

WlUSlUW2US2=WlUW2Ur"1(WlUW2)' 

Since  S^US2  c  r"^(W^UW2),  this  condition  is  also  satisfied. 

Next,  we  have  to  choose  a  nonnegative  integer  function  r  on  £,  that 
is  submodular.  For  W  esd,  we  define  r(W)  =  |w|,  which  clearly  satisfies 
this  requirement  (and  is  in  fact  modular) . 
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We  can  now  apply  Theorem  5.1  to  derive  our  linear  characterization 
of  the  IWS  polyhedron  of  a  bipartite  graph.  To  this  end,  we  set  d  »  1, 
ieV,  in  (5.5),  and  use  the  above  definitions  to  rewrite  (5.6)  as 

(5.6')  x(r"1(W))  +  x(W)  <  |w|,  ¥  W£V2. 

If  we  now  complement  the  variables  x^,  ieV2>  i.e.,  define  new 
variables  x^  *  x±,  itV^  x'  ~  1  -  x±,  ieV2>  then  the  system  (5.5),  (5.6') 
becomes 

(5.7)  0  <  Xi'  <  1,  ieV 

(5.8)  x'(r_1(W))  -  x'(W)  <  0,  vwcv2 

and  Theorem  5.1  asserts  that  the  convex  polytope  P*  defined  by  (5.7),  (5.8) 
has  integer  vertices. 

The  linear  system  of  Theorem  3.1  differs  from  the  above  in  three 
respects.  First,  there  is  an  inequality  (5.8)  for  every  WCVj,  not  just 
those  for  which  W  =  F(S)  for  some  ScV^.  Suppose  that  W  £  T(S)  for  any 
SCVj  and  let  w'  =  r(r_1(W)).  Then  w'cw  and  the  inequality  (5.8)  for  w' 
is  x/(r-1(W))  -  x'(w')  <  0,  which  together  with  (5.7)  implies  the  inequality 

(5.8)  for  W.  Hence  all  such  inequalities  can  be  dropped  without  affecting 
the  integrality  of  the  polytope. 

Second,  (5.8)  does  not  contain  the  inequalities  (3.3)  corresponding 
to  sets  SCV1  such  that  T(S)  =  T(T)  for  some  proper  superset  TcVj  of  S. 

But  if  such  T  exists,  then  the  graph  <  (KjXS)  U  (K^CS))  >  is  disconnected, 
where  <  K  >  is  the  component  of  G  containing  S  and  T,  and  =  KHV^,  i  =  1,2. 
This  is  so  because  the  nodes  in  T\S  /  0  are  not  adjacent  to  any  node  in 
K  \HS)  *  K2\F(T).  As  discussed  in  the  remarks  following  Theorem  3.1, 
the  inequalities  (3.3)  corresponding  to  such  sets  SCV,  are  redundant. 


Third,  the  equation  (3.2)  is  not  present  in  the  system  (5.7),  (5.8). 
This  is  a  genuine  difference  between  the  two  polytopes,  P  defined  by  the 

system  (3.1) -(3. 3),  and  P*  defined  by  (5.7),  (5.8).  However,  the  equation 
(3.2)  defines  a  face  of  P*,  and  since  the  vertices  of  a  face  are  vertices 
of  the  polyhedron,  it  follows  that  P  also  has  integer  vertices.  This  provides 

the  third  proof  of  the  fact  that  P  *  conv  Z. 


6.  Facets  of  the  FMS  Polytope 


In  this  section  we  address  the  question  as  to  which  of  the  inequalities 


defining  the  PMS  polytope  of  a  bipartite  graph  are  essential.  This  is 
obviously  a  matter  of  practical  interest,  as  the  number  of  inequalities  in 
the  system  (3.3)  is  rather  large. 

The  facets  of  a  polyhedron  P  are  its  maximal  (relative  to  inclusion)  non* 
empty  proper  faces.  If  dim  P  Is  the  dimension  of  P,  then  the  dimension  of  a  facet 
of  P  is  dim  P  -  1.  An  inequality  oflc  <  aQ  is  called  facet-inducing  (for  P) , 
if  it  is  satisfied  by  all  x  c  P,  and  the  polyhedron  PH [x|ax  =  aQ}  is  a 
facet  of  P,  i.e.,  has  dimension  dim  P  -  1. 


In  the  remarks  following  Theorem  3.1,  we  have  pointed  out  that  some 


of  the  inequalities  defining  the  PMS  polytope  of  G  are  redundant,  and  that 
the  system  (3.1),  (3.2),  (3.3)  can  in  fact  be  replaced  by  the  smaller  system 
(3.1),  (3.2)  and  (3.3 '),  In  this  section  we  show  that  most  of  the  inequalities 
of  the  latter  system  are  essential,  i.e.,  facet-inducing. 

First,  we  have  to  determine  the  dimension  of  our  polytope.  Let  again 

l  V 1 

P  denote  the  set  of  x  eH1  1  satisfying  (3.1)-(3.3),  shown  in  Theorem  3.1 
to  be  the  IMS  polytope  of  G  »  (V^UVj,  E) . 

The  equality  set  of  the  system  (3.1)- (3.3)  is  the  set  of  those 
members  that  are  satisfied  with  equality  by  all  x  e  P.  A  basis  of  the 
equality  set  is  a  maximal  subset  whose  coefficient  matrix  is  of  full  row 
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For  any  graph  G,  we  define  G,  the  set  of  adjacency  vectors  of  G,  to 
be  the  set  of  all  incidence  vectors  of  pairs  of  nodes  which  are  joined  by 
an  edge.  Thus  G  has  a6  many  elements  as  G  has  edges,  and  each  x  e  S'  has 
exactly  two  components  equal  to  1  and  all  other  components  equal  to  0.  The 
following  Lenina  will  be  useful  in  the  rest  of  this  section. 

Lenina  6,1.  Let  F  be  the  set  of  adjacency  vectors  of  a  forest 
F  *  (V,E)  with  k  components.  Then  F  is  linearly  independent,  | F |  =  )v|  -  k, 
and  every  x  e  F  satisfies  x(K^)  =  x(K^)  for  every  component  (tree)  <  K  > 
of  F,  where  Kj  and  K2  are  the  parts  of  K. 

Proof.  Elementary.  |j 

Theorem  6.2.  Let  X  be  the  set  of  components  of  G  *  (V^UV0,  E),  and 

for  every  <  K  >  eX  ,  let  ■  KOV^,  i  =  1,2.  Then  the  system 

(6.1)  x(Kx)  -  x(K2)  =  0,  *<K>«X, 

is  a  basis  of  the  equality  set  of  (3.1) -(3. 3). 

Proof.  It  is  clear  that  the  equations  (6.1)  are  linearly  independent 

and  belong  to  the  equality  set  of  (3 . 1) - (3. 3) .  Let  F  be  an  edge  maximal 

spanning  forest  of  G,  and  F  the  set  of  its  adjacency  vectors.  Since  every 
pair  of  adjacent  nodes  is  perfectly  matchable,  F  cZ.  By  Lemma  6.1,  F  is 
linearly  independent  and  each  x  e  F  satisfies  (6.1).  Since  |f|  =  |v [  -  k, 
where  k  *  |X|,  no  basis  of  the  equality  set  can  contain  more  than  k  equations. 
But  k  is  the  number  of  equations  in  (6.1),  so  (6.1)  is  a  basis. [j 

Corollary  6.3.  If  G  =  (V^UV2,  E)  has  k  components,  dim  P  *  |v|  -  k. 

Proof.  The  dimension  of  a  polyhedron  inlR^  is  |v|  minus  the  rank 
of  the  equality  set.|| 

We  now  turn  to  the  identification  of  facet  inducing  inequalities. 


The  following  result  will  be  of  use  in  this  task.  We  recall  from  section  2 
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the  definitions  of  7H(G)  as  the  collection  of  perfectly  matchable  node  sets 
of  G,  and  %(£)  as  the  set  of  incidence  vectors  of  such  node  sets. 

Theorem  6.4.  For  any  SCV^,  the  equality 

(6.2)  x(S)  -  x<|-(S)>  =  0 

is  satisfied  by  the  incidence  vectors  of  precisely  those  T  t7!\(G)  such  that 

(6.3)  (V^S,  T(S))  DM  =  0 

for  every  perfect  matching  M  of  <  T  >  . 

Proof.  Let  x  be  the  incidence  vector  of  some  T  cJ^(G).  Clearly,  x 
satisfies  (6.2)  if  and  only  if  |S0t|  ■  |  T  (S)  fl  T  | .  Now  if  (6.3)  holds  for 
at  least  one  perfect  matching  M  of  <  T  >,  then  M  matches  the  nodes  of  SflT 
with  those  of  T(S)  f)T,  hence  x  satisfies  (6.2).  On  the  other  hand,  if  (6.3) 
is  violated  by  some  perfect  matching  M'  of  <  T  >,  then  m'  matches  the  nodes 
of  SflT  with  a  proper  subset  of  the  nodes  of  T(S)  PIT,  hence  <  |P(S)  OtJ 

and  (6.3)  is  violated  by  x.  We  conclude  that  (6.3)  holds  for  at  least  one 
perfect  matching  of  <  T  >  if  and  only  if  it  holds  for  all  perfect  matchings 
of  <  T  >;  and  this  is  the  case  if  and  only  if  the  Incidence  vector  x  of  T 
satisfies  (6. 2). || 

For  any  SCV^,  let  Gg  denote  the  graph  obtained  from  G  by  removing 
the  edge  set  (Vj\S,I"(S)),  i.e.,  let 

gs  -  <  s  ur(s)  >  u  <  (vx\s)  u  (v2\t(s))  >  . 

Then  Theorem  6.4  implies 
Corollary  6.5.  For  any  ScVj, 

%(G)  D  {x|x(S)  -  x(f(S))  -  0}  -  Z( Gs). 
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Theorem  6.4  and  Corollary  6.5  essentially  say  that  for  any  S  EV^, 
the  polyhedron  {x  e  P|x(S)  -  x(r(S))  =  0}  is  itself  a  FMS  polytope,  namely 
the  one  for  the  subgraph  Gg  of  G  obtained  by  deleting  the  edges  in  (Vj\S,r(S)). 

We  are  now  ready  to  state  the  main  result  on  facets  of  conv  2(G) , 
i.e.,  of  P. 

Theorem  6.6.  Let  0  ^  S  cp  V^.  Then  the  inequality 
(6.4)  x(S)  -  x(r(S))  <  0 

is  facet  inducing  if  and  only  if  Gg  has  exactly  one  more  component  than  G. 

Proof.  The  inequality  (6.4)  is  facet  inducing,  i.e.,  the  set 
Pfl  {x|x(S)  -  x(l"(S))  *  0}  is  a  facet  of  P,  if  and  only  if  it  has  dimension 
d  =  dim  P  -  1.  From  Theorem  3.1  and  Corollary  6.5, 

pn  txjx(s)  -  x(r(s))  -  o)  - 
*  conv  2(G)  H  [x[x(S)  -  x(r(S))  *  0}  *  conv  2(Gg). 

From  Corollary  6.3,  dim  P  *  (v|  -  k,  and  dim  conv  2(Gg)  *  jv(  -  kg, 
where  k  and  kg  denote  the  number  of  components  of  G  and  Gg,  respectively. 

Thus  (6.4)  is  facet  inducing  if  and  only  if  kg  ■  k  +  1 . ] j 

At  this  point  there  is  at  least  one  feature  of  Theorem  6.6  that  requires 
immediate  comment.  In  the  remarks  following  Theorem  3.1  we  have  stated  that 
any  inequality  (6.4)  such  that  <  S  U  T(S)  >  is  disconnected,  is  redundant; 
yet  from  Theorem  6.6,  such  an  inequality  may  still  be  facet-inducing,  pro¬ 
vided  that  the  graph  Gg  has  exactly  one  more  component  than  G,  a  condition 
that  is  not  incompatible  with  <  SUF(S)  >  being  disconnected.  So  it  seems 
that  some  facet  inducing  inequalities  are  redundant.  This  is  indeed  the 
case,  due  to  the  fact  that  dim  P  <  |v|,  i.e.,  that  the  equality  set  of 
the  system  (3.1)-(3.3)  is  nonempty.  Every  one  of  the  equalities  satisfied 
by  all  x  c  P  can  be  added,  after  multiplication  with  some  arbitrary  constant. 
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to  any  of  the  inequalities  of  (3.1)-(3.3),  to  yield  another  valid  inequality. 
This  way  infinitely  many  inequalities  may  induce  the  same  facet  of  P,  where¬ 
as  in  any  minimal  linear  system  defining  P,  every  facet  of  P  is  obviously 
represented  by  only  one  (facet  inducing)  inequality.  Thus  we  have  to 
address  the  question  as  to  which  among  the  facet  inducing  inequalities  of 
(3.1)-(3.3)  induce  distinct  facets. 

Before  answering  this  question,  it  will  be  useful  to  restate  Theorem 
6.6  in  the  following  slightly  different  form. 

Theorem  6.6  The  inequality  (6.4),  where  M  S  ^  V^,  is  facet 
inducing  if  and  only  if  G  has  a  unique  component  <  K*  >  such  that  0  t  S*  f  K*, 
where  S*  -  SflK*  and  K*  «  K*^,  i  -  1,2,  and  the  graphs  <  S*UF(S*)  > 
and  <  (K*\S*)  U (K*\T(S*))  >  are  connected. 

This  form  of  the  theorem  (which  can  easily  be  derived  from  the  other 
one)  implies  that  for  all  components  <  K  >  of  G  other  than  <  K*  >,  either 
SDK  i*  0  or  SflK  -  VL. 

Theorem  6.7.  Facet  inducing  inequalities 

(6.5)  x(S)  -x(r(S))  <  0 
and 

(6.6)  x(T)  -  x(T(T))  <  0 

induce  the  same  facet  of  P  if  and  only  if  G  has  a  component  <  K*  >  such  that 

(6.7)  0  *  SflK*  -  TOK*  0  K*HV1. 

Proof .  Since  (6.5)  and  (6.6)  are  facet  inducing,  if  G  has  a  component 
<  K*  >  satisfying  (6.7),  then  K*  is  unique,  and  x  c  P  satisfies  (6.5)  with 
equality  if  and  only  if  it  satisfies  (6.6)  with  equality,  i.e.,  the  two 
inequalities  induce  the  same  facet. 
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Conversely,  if  no  such  <  K*  >  exists,  then  there  exists  ueV^  -S,  veF(S), 
such  that  (u,v)eE  and  either  u,v«TUr(T),  or  u,v*TUr(T).  Then  the  adjacency 
vector  of  {u,v}  satisfies  (6.6)  with  equality,  but  (6.5)  with  strict 
inequality;  i.e.,  the  two  Inequalities  induce  different  facets. || 

We  now  turn  to  the  inequalities  (3.1). 

Theorem  6,8.  The  inequality  xy  >  0  is  facet  inducing  if  and  only 
if  v  is  not  a  cutnode  or  an  isolated  node  of  G. 

Proof ,  If  v  is  neither  a  cutnode  nor  an  isolated  node  of  G,  there 
exists  an  edge-maximal  spanning  forest  F  of  G  in  which  v  has  degree  1. 

Then  the  set  F  of  adjacency  vectors  of  F  contains  a  unique  x  such  that  x^  *  1. 
Therefore,  using  Lemma  6.1,  FU  [0}\(x]  is  a  set  of  dim  P  affinely  independent 
members  of  P,  all  satisfying  xv  *  0.  Thus,  denoting  Q  -  {x  e  p|xv  *  0},  we 
have  dim  Q  >  dim  P  -  1.  On  the  other  hand,  x  c  P\Q,  hence  Q  is  a  proper  face 
of  P;  therefore  dim  Q  •  dim  P  -  1  and  so  *v  >  0  is  facet  inducing. 

Conversely,  if  node  v  is  isolated,  then  x^  *  0  for  every  x  e  P,  and 
thus  x^  >  0  does  not  induce  a  proper  face.  If  v  is  a  cutnode,  let  L  be 
the  node  set  of  a  component  created  by  deleting  v,  and  let  L*  s  LU  {v}. 

Then  every  x  c  P  such  that  xy  *  0  also  satisfies  xfr'flVj)  -  x^'OV^)  *  0. 

But  let  x  be  the  incidence  vector  of  {v,w]  for  any  weL  adjacent  to  v.  Then 
x  satisfies  x(L'HVj)  -  xO/f^)  *  0,  but  x^  +  0.  Hence  the  inequality 
x^  >  0  does  not  induce  a  maximal  proper  face  of  P.|| 

Theorem  6.9.  Facet  inducing  Inequalities  xv  >  0  and  x(S)  -  x(F(S))  <  0 
define  the  same  facet  of  P  if  and  only  if  (vj  ■  K^\S  and  T(v)  cr(K*HS), 
where  K*  is  the  node  set  of  the  unique  component  of  G  satisfying 
0  ft  S  n R*  £  K*  (*  KflVj). 

Proof.  If  the  conditions  hold,  then  the  inequality  xy  >  0  can 
be  obtained  from  x(S)  •  x(T(S))  <  0  by  subtracting  the  equations 
x(Kj)  -  x( Kj)  “  0,  where  •  KflV^,  i  ■  1,2,  for  all  those  components 
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<  K  >  of  G  such  that  K^TiS  t  (9.  Therefore  the  two  inequalities  induce 

the  same  facet.  The  converse  can  be  shown  by  an  argument  analogous  to  the 

one  used  to  prove  the  necessity  of  Theorem  6.7,  and  the  details  are  omitted. |j 

Theorem  6 . 10.  The  inequality  x^  <  1  is  facet  inducing  if  and  only 

if  v  either  has  at  least  two  neighbors,  or  belongs  to  a  two  node  component 

of  G.  In  the  first  case,  no  other  Inequality  (3.1)  or  (3.3)  induces  the 

same  facet.  In  the  second  case,  only  the  inequality  x^  <  1,  where  u  is 

the  other  node  of  the  component  containing  v,  induces  the  same  facet  as 

x  <  1. 
v  — 


Proof . 
define  x  by 


Sufficiency.  If  v  has  two  distinct  neighbors, 

2  if  i  -  v 

,  1  if  i  *  u  or  w 

0  if  i«V\{utv,w}. 


u  and  w. 


Then  xlP,  but  x  satisfies  all  the  constraints  (3.1)-(3.3)  except 
for  the  inequality  xy  <  1.  Therefore  this  Inequality  is  essential,  hence 
facet  inducing,  and  no  other  inequality  of  (3.1)-(3.3)  induces  the  same  facet. 
If  v  belongs  to  a  two  node  component,  with  u  the  other  node,  define 


x  by 

2  if  i  *  u  or  v 

X  ■  i 

1  [o  if  ieV\(u,v}, 

Then  again  x  i  P,  but  x  satisfies  all  the  constraints  of  (3.1)- (3. 3) 

except  for  xy  >  0  and  xy  >  0.  This  shows  that  at  least  one  of  these  two 

Inequalities  is  essential.  But  the  equation  (6.1)  for  the  component  of  G 

containing  u  and  v  gives  xu  -  xy  for  all  x  c  P;  so  x  <  P  satisfies  x^  *  1 

if  and  only  if  it  satisfies  x  ■  1.  Therefore  x  <  1  and  x  <  1  are  both 

v  u  —  v  - 

facet  inducing,  and  they  induce  the  same  facet. 
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Necessity.  If  v  is  an  isolated  node,  x^  =  0  for  all  x  e  P  and  the 

inequality  x^  <  1  does  not  induce  a  nonempty  face  of  P. 

Suppose  now  that  v  has  a  single  neighbor  u,  and  u  has  a  neighbor 
w  ^  v.  If  veV^,  the  inequality  (3.3)  for  S  =  {v}  is  xv  -  xu  <  0,  or  xv  <  xu- 

If  veV2»  this  same  inequality,  though  not  part  of  (3.3),  can  be  derived  as 

the  inequality  (3.13)  for  T  *  {v}.  Therefore  every  x  e  P  that  satisfies 
x  *  1  also  satisfies  x  *  1.  But  the  converse  is  not  true,  since  the 
adjacency  vector  x  of  [u,w]  belongs  to  P,  while  x^  “  1,  x  *  0.  Therefore 
the  inequality  xv  <  1  does  not  induce  a  maximal  proper  face  of  P.|| 

From  the  last  four  theorems  it  follows  that  the  set  of  constraints 
(3.1),  (3.2),  (3.3 ')  comes  very  close  to,  though  is  generally  not  exactly, 
a  minimal  linear  system  defining  P.  To  make  it  minimal,  one  has  to  remove 

-  every  inequality  xv  >  0  such  that  v  is  either  an  isolated  node 
or  a  cutnode; 

-  every  inequality  <  1  such  that  v  has  less  than  two  neighbors  and 
does  not  belong  to  a  two  node  component;  and,  finally, 

-  every  Inequality  x(S)  -  x(r(S>)  <  0  such  that  |K*\S|  *  1  and 
T(K*\S)  C  T(K*  H S) ,  where  K*  is  the  unique  component  of  G  such 
that  0  i  K*0S  i  K*  (-  K*nvx). 

This  still  leaves  a  large  number  of  inequalities,  that  can  be 
exponential  in  the  size  of  G.  The  following  example  illustrates  the 
contents  of  this  section  and  also  is  a  case  where  the  minimal  defining 
system  for  P  is  exponential. 

Let  Gq  be  the  graph  of  Fig,  6.1,  consisting  of  n  pairs  of  nodes 
£ut,  v1},  each  pair  joined  by  an  edge,  plus  a  node  uQ  adjacent  to  every 
v^,  i  ■  l,...,n,  and  a  node  vQ  adjacent  to  every  u^,  i  «  l,...,n.  Let 
V1  *  &V  . uQ},  V2  -  {vo,  v1#...,vnJ. 
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Using  the  results  of  this  section,  we  obtain  the  following  minimal 
defining  system  for  P: 


0  <  x  <  1,  u«V, 

—  u  —  ’  1 

(6.14) 

0  <  x  <1,  vcV„ 

—  v  —  *  2 

(6.15) 

x(Vj)  -  x(V2)  =  0 

(6.16) 

X(S)  -  x(r(S))  <  0,  V  S:0  +  S 

(6.17) 

xu  *  x<v2Xtvo})  ^  °' 

0 

The  Inequalities  (6.14)  and  the  equation  (6.15)  are  easily  seen  to 

be  needed.  For  any  nonempty  S c {»,... ,u  },  v  *r(S),  so  <  S  U  F(S)  >  Is 

j.  no 

connected;  and  uo«V^\S,  so  <  (V^\S)  U (V2\T(S))  >  is  also  connected;  hence 
from  Theorems  6.6  and  6.7,  the  inequalities  (6.16)  all  Induce  distinct 
facets  of  P.  Further,  since  S  ^  Vj\[uo)  implies  jVj\s|  >  2,  the  facets 
induced  by  these  inequalities  are  also  distinct  from  those  induced  by  any 
of  the  inequalities  (6.14)  (Theorem  6.9).  Finally,  since  <  [uQ} U (V2\{vq})  > 
and  <  (Vj\{uq})  U  {vQ}  >  are  connected,  inequality  (6.17)  defines  a  facet  of 
P,  which  is  easily  seen  to  be  distinct  from  the  facets  induced  by  any  of 
the  other  inequalities. 
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It  remains  to  be  shown  that  the  omission  of  the  reuu-.aing  inequalities 

of  (3.3)  is  justified.  If  S  *  Vj\{u  },  then  from  Theorem  6.9  the  inequality 

(3.3)  induces  the  same  facet  as  x  >0.  Now  let  u  eS. 

u  —  o 

o 

If  S  *  {uq},  we  tow  the  inequality  (6.17).  Now  let  S  /  Cu03-  Then 
r(S)  »  V^.  so  <  (Vj\S)  U  (V^XIXS) )  >  is  connected  if  and  only  if  |Vj\Sj  <  1. 

If  \V1\S\  ■  0,  then  S  *  and  the  inequality  (3.3)  is  Implied  by  the  equa¬ 
tion  (6.15).  If  Vj\5  *  {u^}  for  some  ie{l,...,n},  then  from  Theorem  6.9 

the  inequality  (3.3)  induces  the  same  facet  as  x  >0.  This  covers  all 

ui 

the  cases. 

Notice  that  the  number  of  inequalities  (6.16)  for  G  is  2°  -  2,  hence 

n 

exponential  In  n. 

Although  the  number  of  inequalities  in  our  linear  characterization  of 
the  PMS  polytope  of  a  bipartite  graph  may  be  large,  this  characterization  is 
still  computationally  useful.  Indeed,  a  linear  program  whose  constraint 
set  includes  the  system  (3.1)-(3.37)  can  be  solved  by  generating  the  in¬ 
equalities  (3.3')  as  needed.  However,  the  development  of  such  a  procedure 
goes  beyond  the  scope  of  this  paper. 
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sets  of  G.  We  derive  this  result  by  three  different  approaches,  using  linear 
programming  duality,  projection,  and  lattice  polyhedra,  respectively.  The 
projection  approach  is  used  here  for  the  first  time  as  a  proof  method  in  poly¬ 
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inequalities  of  our  linear  defining  system  from  the  redundant  ones. 


